Growth enhanced surface diffusion and elastic instabihty on amorphous solids 
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A continuum model for growth of solids is developed, considering adatom deposition, surface 
diffusion, and configuration dependent incorporation rate. For amorphous solids it is related to 
surface energy densities. The high adatom density leads to growth enhanced dynamics of (a) Mullins' 
classical equation [J. Appl. Phys. 28, 333 (1957)] without, and (b) of the Asaro-Tiller-Grinfeld- 
Srolovitz instability with lateral stress in the growing film. The latter mechanism is attributed to 
morphologies found in recent experiments. 
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The theoretical approach to kinetic roughening [Q, 
due to its roots in Statistical Mechanics, relates vari- 
ous growth processes on solids to so called universal- 
ity classes. Their distinction allows to identify essen- 
tial properties for the large scale morphology beyond the 
microscopic mechanisms, say certain symmetries or con- 
servation laws, which manifest themselves in terms of 
universal scaling exponents. 

One such important property is mass conservation sup- 
pressing adatom loss to the surrounding space, which al- 
lows for surface relaxation only through transport along 
the surface. The total solid mass then is always equal to 
its initial value plus the integrated incoming flux. 

Molecular Beam Epitaxy (MBE) garantees such condi- 
tions, and growth parameters like for instance flux inten- 
sity, temperature, chemical composition of the adsorbate 
can be easily and precisely controlled It has been 
widely used and refined for growth of semiconductor and 
metal crystals. However, under typical conditions the dif- 
fusion length and average crystal step distance are quite 
large, imposing their own features to the morphology 
such that any theoretically predicted asymptotic 
scaling regime lies beyond the reach of experimental ob- 
servation. 

On the other hand it may become observable in grow- 
ing amorphous substances where it seems plausible that 
intrinsic lengthscales remain small. Some experiments 
have been performed in recent years looking for kinetic 
roughening on amorphous solids ||5|-|l0t| . It turns out that 
their theoretical interpretation in terms of standard con- 
tinuum [0 and discrete ||l^,|l3| models leaves two main 
open questions which are addressed in this letter. First, 
far from equilibrium the kinetics on the surface are dif- 
ferent from relaxation by thermally activated adatoms 
leading to Mullins' continuum equation ]ll| . Second, in 
some experiments there appears a very pronounced in- 
termediate lateral lengthscale together with a mounding 
instability (^-|^ whose origin was unclear. Obviously it 
cannot be due to step edge barriers as is the case on crys- 
tal surfaces [plH . Another suggested mechanism, deflex- 



ion of incoming particle trajectories or "steering" Q, is 
in some cases relevant at grazing incidence but negligible 
for normal beams [ p^ . 

Here a new continuum approach to growth with surface 
diffusion is developed resulting in a flux dependent coef- 
ficient for surface evolution (see Eqs. (H) and (Jl2|)), the 
counterpart to equilibrium relaxation Eq. (^ pT] . The 
mounding instability is shown to be due to lateral elastic 
stress in the growing solid |5|P, p^ , p^ , which is respon- 
sible for the well known Asaro-Tiller-Grinfeld-Srolovitz 
(ATGS) instability ||l|-|2|]. Surface modulations allow 
for relaxation of a laterally stressed solid. For large 
enough wavelengths this energy gain overcomes the cost 
in additional capillary energy, disfavouring and destabi- 
lizing a flat surface. Due to the nonequilibrium kinetics 
it appears in a new light: Its driving force is an equi- 
librium, the dynamical evolution a nonequilibrium phe- 
nomenon. In the following calculations the surface ener- 
getics account well for the observed wavelength (the size 
of emerging mounds), and the nonequilibrium kinetics 
yield the linear growth rate of the instabilty. 

To begin, the continuum model of equilibrium relax- 
ation is briefly recalled [0. Surface configurations are 
described by a space and time dependent height field 
i/(x, i), neglecting overhangs and voids. The driving 
force for equilibrium relaxation is the surface free energy, 
J- — J (Pxj^yi + {VHy. If an atom is added somewhere 
at the surface this total area may change, resulting in a 
extra energetic contribution /i = —VL'yV^H ||2^]. /i is 
often called chemical potential and Vl = denotes the 
atomic volume. Spatial variations in bias the adatom 
diffusion and create a macroscopic mass exchange. Lo- 
cal mass (and volume) conservation impose a continuity 
equation for the surface dynamics, which for <C ksT 
Ip3| reads 



dtH 



(1) 



A reasonable assumption for the adatom diffusion coef- 
ficient is D = a^/ro 6"^°'''/'^^^, an Arrhenius form with 
some activation barrier ScS (which on a random surface 
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comes from an average over waiting times at each site 
iQ) and an "attempt frequency" I/tq Eq. (P pre- 
dicts an exponential decay with rate —nk^ for fluctua- 
tions with wavenumber k, which has been veryfied ex- 
perimentaUy On a phenomenological basis Eq. (|l|) 
has also been applied to growth with relaxation by cur- 
vature driven surface diffusion, where under deposition 
("shot") noise it results in a power spectrum of height 
fluctuations ^(k,*) = {\H(k,t)\^) oc g. Of course 

under growth the coefficient k takes a different form as 
above in Eq. (|l|) which is adressed next. 

On a growing solid there is a constant supply of mobile 
adatoms out of the beam hitting the surface. The number 
of thermally activated adatoms, responsible for equilib- 
rium relaxation, now becomes negligible. Here an ansatz 
is presented where the incorporation rate I{H} depends 
on the surface configuration i?(x, t) and connects it thus 
to the adatom density p{x, t), 



DF 



dtH = I{H}np 
dtp^V ■D{H,p}Vp~I{H]p 



(2) 



Note that H + flp is a conserved quantity, increasing by 
the average growth velocity F. The problem is to find 
the dependence I{H}, but before turning to that point 
its role in a phenomenological equation as (|l|) shall be 
worked out. 

Eqs. (|^) are expanded around the "fiat" growing sur- 
face Ho{t) = Hq + Ft, on which the incorporation rate 
takes some value Iq , and around the average adatom den- 
sity po = F/{QIo). Small deviations /i(x, t) = H{x,t) — 
Ho{t) and c(x, i) = ^2(p(x, t) — po) obey to hnear order 



F 

dth ~ — — /i *h + loc 
lo 

dtc =—h*h+ {D\/^ - Io)c, 



(3) 



where D = D{Hf){t), po}. I{H} — Io ~ Ii * h is approx- 
imated linearly (not necessarily local in space) and will 
be evaluated in its Fourier transform Ii (k) . 

In the long wavelength limit, Dk^ <C Iq (beyond the 
diffusion length £d, the typical distance an adatom moves 
before incorporation), and also |/i(k)|i^//o ^ Iq the 
eigenvalues of (0) are 



n F 

A2 = —Iq 



(4) 



with relative errors ©(DP/Jo) + 0(i^|4 (k)|//2). A2, re- 
maining finite in the long wavelength limit, rules the en- 
hanced (diminished) incorporation under higher (lower) 
adatom density. The eigenvector of Ai mainly lies in di- 
rection of h, so in the considered large scale limit the 
linearization of Eqs. (Q) turns into 



dtHKt) = --^k'hik)h{k,t), 



(5) 



a non-equilibrium version of Eq. 

On crystal surfaces Ii(k) depends essentially on the 
step configurations, but on amorphous solids it is pro- 
portional to variations in energy density at the surface, 
and can be derived in a mean field type approach. For a 
diffusing adatom the amorphous surface consists of po- 
tential wells of various depths. They have an average 
distance a from each other and a probability distribu- 
tion of depths n{e/eo)/£Q, where the energy scale Eq is 
explicitely included in the notation. In a simple view 
of the nonequilibrium growth process an atom gets de- 
posited, diffuses, and is finally incorporated into the solid 
in a sufficiently deep potential well. Dimer formation 
or other types of nucleation can be neglected if sticky 
sites are denser than adatoms — a condition to be ver- 
ified a posteriori. An adatom "sticks" to a site, if it 
cannot escape until it is buried by the further growing 
solid. Thus the depths e of "sticky" potential wells fulfill 
tqc'^/^'^'^ > a/F. Assuming an exponential distribution 



of energy depths, n(e/eo) = exp(— e/eo) |28|, this yields 
a relative proportion of sticky sites 



'-d 



ks T/ eo 



(6) 



which also defines the diffusion length i^, because the 
fractal dimension of a random walk is two and the adatom 
can "fully explore" its surroundings on the surface. The 
average incorporation rate is the inverse of the time 
needed to hit a sticky site, so /q = D/£^. 

As seen in the context of Eq. (|l|) deviations from a flat 
surface change the energy density and the chemical po- 
tential near the surface. A natural way to account for 
this is shifting the energy scale Eq Eq — p ^ Eq — fl£. A 
positive additional energy density £ lowers the depth of 
potential wells encountered by the adatoms, and changes 
the incorporation rate by 



a kBTnS 
I -h = -Iq log 

Ftq Eq Eq 



^2/0 log . (7) 

a Eq 



The change is proportional to the relative energy change 
for adatoms. Flux and temperature enter via id and /q. 

Now the surface energetics have to be evaluated in or- 
der to obtain /i(k) in Eq. (|^) via (^. As above in Eq. 
(|l|) also here the change in surface free energy enters, so 
there must be a capillary contribution |23] 



£^i"(x,i)--7V2Mx,<). 



(8) 



The other important part of £ comes from elastic de- 
formations of the growing solid. In the experiments con- 
sidered here amorphous metallic glasses (i.e. the alloy 
Zrg5Al7.5Cu27.5) grow under lateral expansive stress the 
origin of which is not fully clarified It builds 
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up during growth and reaches a constant level at a film 
thickness of about 50 nm. The order of magnitude is 
then 1 GPa, corresponding to a lateral deformation of 
a Ri 0.1 to 1 %, given Young's modulus to be roughly 
E w 10^ GPa Applying the linear relations of 

stress and strain for small deformations one can cal- 
culate the strain and stress tensors in the film with a flat 
surface, and a^j, as well as the elastic energy density 

Changes in shape are slow compared to mechanical 
balancing inside the body, so the strain and stress fields 
follow the surface configuration quasistatically. Surface 
variations can be seen as perturbing (t°- by an addi- 
tional stress field Tij with boundary conditions r^^ — 
—dihEa/ {1 — a) for i — x,y and Tzz = to linear order 
in V/i (cr without indices denotes the Poisson number). 
Green's function corresponding to the geometry of the 
body |29j yields the perturbative and corresponding 
strain Wij throughout the solid. In view of Eq. only 
the corrections to the elastic energy density at the 
surface to linear order in nj and Wij (and therefore V/i) 
are needed, 



a 



E 



TT 1 - 2cr 



Curved parts of the surface also compress or elongate 
the solid, contributing to the energy density by 



1-2(7 



aj V^/i(x, t). 



(10) 



The derivation of these elastic energy densities will be 
presented elsewhere in more detail [ 0O|. Here the next 
step is to gather Eqs. (||), (||), and ( |10D to the Fourier 
transform of the energy density variations at the surface. 



£'"\k,t) = B{k) h{k,t) 



(11) 



1-2^ \ ,2 E 2l+^|.| 

aa— 7fc a |k| 



1-a' 



h(k,t), 



which is local in k. With i?(k) defined above, the lin- 
ear surface evolution in Eq. (^ including capillary and 
elastic effects becomes 



dthik,t) 



-2£2log^F^^ffifc2;^(k,i). 

a So 



(12) 



In particular, without any elastic effects the kinetic co- 
efficient in Eq. (0) is K = 2£l\og{£d/a)F{n'y/eo). The 
full -B(k) reflects the ATGS instability: It is positive 
for large |k|, but negative below a critical wavenum- 
ber kc, attaining a minimum at the wavevectors |k| = 
A;* =3/4a2/7r(l-|-cr)/(l-CT)£'/7, where B{k) takes the 
value —"fkl/3. In Eq. because of (even small) ini- 
tial roughness and noise in the deposition and diffusion 
processes from this linear instability random patterns 



of buckles with typical scale A* = 27r/fc* emerge grow- 
ing exponentially in amplitude [ pl| with a rate 1/t* = 
2/3 £l \og(ej a)F{n-f/eo)kt = nkfjt. 

So far a nonequilibirum dynamic equation for grow- 
ing surfaces, Eq. (|l^), has been derived perturbatively 
close to a horizontal interface. Now it is compared to 
experimental results p|-f7|,p|,p^ . The instability has been 
carefully observed e.g. for Zr65Al7.5Cu27.5 but the 
lateral stress is documented only for related materials 
1^,^. The best test would be to measure both simulta- 
neously, but already with the given information an order 
of magnitude estimate for A* and 1/t* is possible. 

Here are the experimental parameters: The elastic con- 
stants were given above, the surface tension is 7^2 J/m^ 

f , temperature fcBT«4xlO^^^ J, growth velocity F — % 
/s, atomic size a = Q}/'^ Ki 3 A Order of magni- 
tude estimates are used for the time between adatom 
hop attempts, tq ~ 10^^^ s and surface energetics, 
eo~10^^^ J (about 1 eV). This yields a diffusion length 
td~^ nm, which suits well the observation that atoms 
move a few diameters before incorporation Indeed 
this is smaller than the average distance between mobile 
adatoms y^Q/po ~ £d\/D/F w a^Jso/kBT (D is calcu- 
lated from the distribution of waiting times in "unsticky" 
sites [|4|), so dimer formation is suppressed. Given these 
values the theoretical predictions nm and 

l/nwlO-2 1/s. 
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m thickness (nm) 



(dashed lines) to experi- 



100 200 

film thicl<ness (nm) 

FIG. 1. Comparing 1/r, and A» ^ 
ments Q (diamonds). Left: mean square roughness, grow- 
ing exponentially after thickness 30 nm. Dotted line shows 
growth rate 6.5xl0~'^ 1/s. Right: lateral buckle size. 



This fits fairly well to experimental observations, as il- 
lustrated in Fig. 0. The growing film develops buckles, 
which after a thickness of about 30 nm take a constant 
lateral size of i?c ~ 17 nm. From 30 to 240 nm film thick- 
ness their vertical amplitude increases exponentially with 
rate 6.5 xlO""^ 1/s. Surprisingly this quantitative expres- 
sion of a linear growth instability is not adressed explic- 
itly in the original work where the authors focus on 
an observed early time algebraic increase for both quanti- 
ties. It is caused by kinetic roughening of large k modes 
at times before r, pc| ]. Besides, in experiments on re- 
lated materials the lateral homogeneous strain a reaches 
a constant level only after about 50 nm film thickness 
[p|, p^. S o only then the linear instability as described by 
Eq. (|12|) with Bik) constant in time becomes visible. 

Some remarks comparing to different interpretations 
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of the observed instability are in order. First, the 
elastic energy density in Eq. ( 9) is calculated for film and 
substrate of the same material ||2^,0 . This does apply to 
the experimental system, where the Zr alloy film was de- 
posited on a previous 100 nm thick layer of the same ma- 
terial. In particular, effects of perfectly rigid substrates 
p2[ won't be observed. Second, the partial relaxation 
of the film close to the modulated surface will not pro- 
duce a measurable relief of total stress in the layer. Even 
close to the surface the stress is lowered only by a factor 
l—0{\Vh\) |30j and the method of substrate deformation 
measures only an average stress across the whole film . 
Third, compression of convex expansion of concave parts 
by surface tension should not be measurable in the to- 
tal stress changes, since positive and negative curvature 
compensate each other. Besides it is only a minor effect 
(compare Eq. (^0|) to (||) and (^). Forth, recently a local 
continuum equation with linear terms dth = {vk^ ^ Kk'^)h 
has been fitted to experimental results 0. ft would be 
interesting to see whether a destabilizing term oc |k|^ as 
in Eq. (||) can give a better description. 

In conclusion, in this letter a theoretical framework for 
a continuum thcry of surface growth with diffusion has 
been constructed. As in the fundamental lattice mod- 
els [|l^,|3| the basic processes are particle deposition and 
diffusion until an energetically favourable site is reached. 
Surface free energy stabilizes the interface by a configu- 
ration dependent attachment rate. For amorphous solids 
a mean field type of approach yields Eq. a configu- 
ration dependence through the energy density near the 
surface, resulting in a nonequilibrium evolution equa- 
tion (p^. Second, an experimentally observed growth 
instability on amorphous films has been shown to be 
stress induced. Its spatial properties can be explained by 
standard energy arguments of the Asaro- Tiller- Grinfeld- 
Srolovitz instability, its temporal evolution needs the 
above nonequilibrium framework. Far from equilibrium 
the instability is enhanced. A challenge is the extension 
to crystal growth, where steps act as sticky sites and can- 
not be treated in a simple mean field way, and where is- 
land nucleation becomes important. Random nucleation 
leads to different step configurations on maxima com- 
pared to minima, so it contributes to Ii (k) which may be 
a way to understand the relation noci'^F obtained from 
dimensional analysis |2I 3^, where £d denotes the av- 
erage distance between island nuclei on crystal terraces, 
conceptually different from £d here. Also the elastic in- 
teractions should be worked out for crystals and for het- 
eroepitaxy with different rigidities of substrate and grow- 
ing film, which would enable very important applications 
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